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Thermosonde  and  UHF  Radar  Measurements  of  Cn 
at  Westford,  Massachusettes- July  1981 


1.  INTRODUCTION 

In  July  1981  over  30  specialized  AFGL  radiosondes  were  launched  from  the 

Millstone  Hill  radar  facility  in  Westford,  Massachusetts.  1  The  purpose  was  to 

provide  turbulence  data  to  determine  the  statistics  of  refractive  index  fluctuations 

and  their  influence  on  optical  propagation  degradation.  The  published  data  and 

models  are  insufficient  to  predict  the  role  of  turbulence  and  the  extent  turbulence 

degrades  the  performance  of  new  communications,  guidance  and  navigation  systems. 

The  particular  objectives  of  these  experiments  were:  (1)  to  expand  the  limited  data 

2 

base;  (2)  to  compare  the  altitude  profiles  of  C  from  measurements  by  thermosonde, 

n  2 

radar,  aircraft,  and  stellar  scintillometer;  (3)  to  estimate  an  average  profile 
for  day  and  night  conditions;  and  (4)  to  estimate  coherence  scales  based  on  the 
thermosonde  measurements. 

(Received  for  publication  22  February  1984) 

1.  Brown,  J .11.,  Good,  R.  E. ,  Bench,  P.  M. ,  and  Faucher,  G.  E.  ( 1982)  Sonde 
Experiments  for  Comparative  Measurements  of  Optical  Turbulence, 

At'GL  - T R- 8 2  - 007 9,  AD  A1  18740. 
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2.  PROGRAM  DESCRIPTION 


Coordinated  experiments  were  conducted  near  the  Haystack  observatory  at 
Millstone  Hill,  Westford,  Massachusetts  from  10  February  198 1  until  8  August  1 981 .  1 
Turbulence  was  measured  by  Al'Gl.  modified  radiosondes  and  thermosondes,  the 
Northeast  Kadio  Observatory  Corporation  (Nh’ROC)  440  MHz  UHF  scatter  radar, 
an  AFWL  stellar  scintillometer,  and  a  contractor  instrumented  aircraft.  Per¬ 
sonnel  from  the  Northeast  Radio  Observatory  Corporation  operated  the  radar.  The 
airctaft  mounted  micro-thermal  probe  measurements  were  provided  by  Airborne 
Research  Associates,  Bedford.  Massachusetts.  AFG L  personnel  launched  the 
radiosonde /thermosonde  combination  and  conducted  the  scintillometer  measurements. 

The  plan  was  to  operate  all  instruments  simultaneously  but  this  was  not  possible. 
Personnel,  scheduling,  and  weather  conflicts,  to  a  great  extent,  dictated  when  each 
instrument  was  available.  Consequently,  experiments  were  conducted  with  two  or 
three  of  the  instruments:  thermosondes  and  radar;  radar  and  scintillometer; 
thermosondes,  radar,  and  aircraft. 


3.  THEORY 


3.1  Thermosonde 


The  AFGl.  thermosonde1  is  an  ac  Wheatstone  bridge  that  measures  small  resis¬ 
tance  differences  between  tw'o  unheated  tungsten  wires  exposed  to  the  ambient 
atmosphere.  A  2-m  weather  balloon  carries  the  instrument  to  altitudes  up  to  30  km, 
and  the  downlink  transmission  is  facilitated  by  a  modified  VIZ  Corp.  radiosonde. 

The  probe  wires  are  separated  horizontally  by  1  m  and  the  difference  in  resistance 
is  translated  to  a  difference  in  temperature  by  laboratory  calibration  and  measure¬ 
ments  of  the  temperature  coefficient  of  resistance  a  of  the  tungsten  probes.  Specific¬ 
ally,  R  =  R  +  a  R  (T-  T  ),  where  a  is  evaluated  at  temperature  T  =  0°  C.  The 
■’pooo  2  V  8 

inertial  range  refractive  index  structure  parameter,  is  calculated12  from  the 

expression: 


C"  ( z > 
n 


79.9  X  1C 


P(z> 
~~2 - 

T  (z> 


2 

C*  (z) 


where  P(z)  is  pressure  in  mbar  and  z  is  altitude.  CTp  is  the  temperature  structure 
parameter,  and  it  is  evaluated  from  the  mean  square  temperature  fluctuation  as 
provided  by  the  thermosonde.  Pressure,  temperature,  humidity,  and  CVp  are 


2.  Tatarski,  V.  [.  11961)  Wave  Propagation  in  a  Turbulent  Medium,  McGraw-Hill  Co 
New  York. 
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sampled  by  the  radiosonde  and  transmitted  once  every  2  sec  or  about  every  10  ni 
altitude,  \ltitudo  could  be  obtained  independently  from  the  tracking  data. 

3.2  Radar 

The  Millstone  Hill  440-IUHz  radar  was  operated  bv  Northeast  Radio  Observatory 

3  2 

Corporation  personnel  to  measure  continuous  altitude  profiles  of  C  .  fluctuations 

in  the  atmospheric  index  of  refraction  arising  from  gradients  in  temperature  and 
humidity  generated  by  turbulence  were  the  source  of  signal  returns  from  the  sensi¬ 
tive  l  111  radar.  Turbulence  is  evaluated  from  the  wind  shifted  Doppler  frequency 
spectrum  of  the  signal. 

2  4 

The  radar  volume  reflectivity  71  is  related  to  C  bv 

1  n  • 

■>  j  3 

G"  -  II)  10.  331  3 
n 

where  >  is  the  radar  wavelength.  The  scattering  is  assumed  to  arise  from  homo¬ 
geneous  isotropic  turbulence  that  fills  the  entire  scattering  volume.  The  reflectivity 
is  related  to  the  measured  signal  to  noise  ratio,  SNR,  by ^ 

i)  =  R^  ^  (SNR) 

where  R  is  the  radar  calibration  constant,  r  is  the  radar  slant  range,  and  Ar  is 

c 

the  range  gate  resolution.  The  radar  azimuth  and  elevation  angles  were  aligned  as 
clo.-.e  as  possible  with  the  radiosonde  path,  but  not  so  close  as  to  cause  interference. 
In  these  experiments,  the  radar  had  a  height  resolution  of  about  1  km  with  each 
profile  being  evaluated  approximately  every  minute.  Data  below  5  km  were  eliminated 
so  as  to  eliminate  large  ground  clutter  noise  and  to  avoid  high  sensitivity  to  water 
vapor.  Data  above  18  km  were  not  reliable  due  to  the  small  SNR  in  this  region. 

3.3  Stellar  Scintillometer 

The  stellar  scintillometer  is  a  36 -cm  Schmidt-Cassegrain  telescope  that  focuses 
starlight  into  an  attached  instrument  assembly.  The  scintillation  intensity  is 
measured  at  three  different  spatial  wavelengths  between  5  and  15  cm.  These  spatial 


3  Hood,  l{.  E.  ,  W  atkins,  B.  J.  ,  Quesada,  A.  I-'.  ,  Brown,  J.  If.  ,  and  Lo riot,  G.  B. 
i  1982)  Radar  and  optical  measurements  of  C^.Appl.  Opt.  ,  2U  3373-3376. 

4.  Ottersten,  11.11.  (1969)  Radar  backscattering  from  the  turbulent  clear  atmos¬ 

phere,  Radio  Science,  4:125  1-  1255. 

-  VSA, 

5.  Balslvv,  B.  B.  ,  and  Gage,  K.  S.  (1980)  The  MS T  radar  technique  potential  for 

middle  atmospheric  studies,  Bure  Appl.  Geophvs.  ,  118:452. 


wavelengths  are  linearly  combined  to  form  seven  composite  path  -  weighting 
functions  to  provide  seven  different  altitude  regions  of  the  turbulence  measure¬ 
ments.  Halfwidths  of  the  weighting  functions  increase  from  1.4  km  at  2.25  km 

2 

altitude  to  6.  0  km  at  18  km  altitude.  C  is  calculated  from  the  normalized  stellar 

n 

irradiance  variance  detected  through  the  spatial  filter, 

1/2 


2  2 

o"  where  o.  ~ 

If  lf 


[In  (I.  +  1,1  -  In  (I„  +  I 


I 


and  1.  and  I  are  the  mean  and  fluctuating  part  of  the  signal  passed  by  the  filter 

I  1  j 

reticle  and  T7  and  1^  are  the  mean  and  fluctuating  part  of  the  signal  reflected  by  the 


filter  reticle. 


3.4  Aircraft  Measurements 

Airborne  Research  Associates  (ARA),  under  contract  to  Al-'GL,  equipped  an 

aircraft  with  4.5  /am  diameter  heated  tungsten  wires.  These  were  mounted  on  a 

wingtip,  one  above  the  other  with  a  separation  of  0.85  m.  Wheatstone  bridge 

measurements  of  the  voltage  difference  between  the  two  wires  were  recorded  along 

with  simultaneous  measurements'  of  air  temperature,  pressure,  humidity,  and 
2 

airspeed.  C'T  was  calculated  from: 


C2  =  |  (AT  >2 
T  I  rms  t 

where  N  is  the  noise  level  with  contributions  from  the  bridge  and  from  velocity  fluc¬ 
tuations.  The  noise  level  was  estimated  by  measuring  the  bridge  voltage  under  rela¬ 
tively  quiet  signal  conditions.  Additionally,  the  contribution  due  to  velocity  fluctua¬ 
tions  was  sensed  by  a  single  hot  wire  and  a  constant  temperature  anemometer  bridge. 
2*2 

C  was  then  calculated  from  C'  in  the  same  manner  as  for  the  thermosonde, 
n  T 


4.  RAWINSONDE  RESULTS 


4.1  Temperature  and  Humidity 


Ambient  temperature  was  measured  by  a  thermistor  rod  mounted  on  a  plastic 
outrigger  attached  to  the  radiosonde.  Relative  humidity  was  measured  by  a  carbon 


6.  Ochs,  G.  R.  ,  Wang,  T.  ,  L.awrence,  R.  S.  ,  Clifford,  S.  F.  (1976'  Refractive- 

turbulence  profiles  measured  by  one-dimensional  spatial  filtering  of  scin¬ 
tillations,  _Aj3£L_0£t.  ,  ^5^2564-2510. 

7.  i  airall,  C.  W. ,  and  lilarkson,  R.  (1979)  Aircraft  Measurements  of  Micro- 

Meteorological  Parameters  at  Panama  City,  Florida  in  1978.  Naval 
Postgraduate  School  Technical  report. 
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hygristor  mounted  in  a  flow  duct  inside  the  radiosonde.  These  measurements  were 
sampled  every  2  sec  or  about  every  10  m.  Figures  1-19  are  profiles  of  temperature 
and  relative  humidity  for  19  flights.  Since  the  hygristor  is  insensitive  to  relative 
humidity  below  20  percent,  the  relative  humidity  plots  are  cut  off  at  10  km  where 
this  condition  normally  prevails. 

4.2  Wind  Speed  and  Direction 

A  technique  based  upon  radio  signal  phase  shift  was  used  to  measure  the  balloon's 
position  as  a  function  of  time.  Derivatives  of  the  position  are  evaluated  to  estimate 
wind  speed  and  direction.  Figures  20-36  are  profiles  of  the  wind  speed  and  direc¬ 
tion  of  17  flights.  Outlier  points  were  removed  and  the  raw  data  smoothed  and 
filtered  by  a  low  pass  filter  whose  cutoff  was  about  0.  03  Hz. 

Mean  square  wind  speeds,  W^,  are  evaluated  for  these  flights  for  use  in  the 
8  2  2" 

Hufnagel  model.  Values  of  W  are  listed  in  Table  1. 

4.3  Thermosonde  Cn2  Profiles  -  Hufnagel  Comparisons 

2 

Weighted  500-m  log-averaged  of  the  10  m  samples  of  thermosonde  C  data  are 
plotted  in  Figures  37-50.  The  weighted  average  is  calculated  from 


,.2  .  .  n  .  - 1 

Cn  (z)  =  ^  log 


i  y* 


where  N  is  the  total  number  of  points  over  Az  =  500  m,  and  n  is  the  number  of  points 

above  the  thermosonde  noise  limit  of  0.  002°  C.  The  noise  limit  is  reached  in  regions 

2 

of  very  low  turbulence.  A  logarithmic  average  is  performed  since  C'n  is  represented 

by  a  log-normal  distribution.  The  mean,  multiplied  by  or  divided  by  one  standard 

deviation,  is  shown  on  either  side  of  the  mean.  The  curves  are  broken  and  replaced 

by  circles  where  radio  loss  of  signal  occurs  and/or  where  the  number  of  data  points, 

2 

n,  were  small.  Actual  instantaneous  C  data  are  plotted  for  the  first  500  m  above 

n  2 

the  surface  since  in  this  boundary  layer  region,  C  falls  off  rapidly  with  altitude  and 
a  500  m  average  is  not  representative.  Obviously,  the  profiles  are  space-  and  time- 
dependent;  however,  the  trend  curve  is  considered  valid  over  the  length  of  the 
experiment. 

2  fi 

To  estimate  the  C  profile  from  rawinsonde  meteorological  data,  Hufnagel 

developed  a  mean  square  wind  speed  dependent  model.  This  simple  empirical  model 

depends  only  on  W  ,  the  mean  square  wind  speed  over  5  km  to  20  km.  For  altitudes 

of  3  km  to  24  km  the  model  is  given  as, 

8.  Hufnagel,  K.  F.  (1974)  V  ariations  of  atmospheric  turbulence,  The  Infrared 
Handbook,  liSGPO,  Washington,  D.  C.  ,  Chapter  6,  pp  1-56. 
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\  alues  of  W  are  listed  in  Table  1.  In  three  eases,  QX45,  QX35,  and  QX43,  wind 

data  are  not  obtained  up  to  20  km;  nevertheless,  \\  calculated  to  the  maximum 

altitude  is  presented  for  comparison.  Profiles  calculated  with  the  Hufnagel  model 

are  depicted  in  Figures  37-50  as  solid  dark  lines.  In  most  cases,  order-of- 

2 

magnitude  agreement  is  found  between  the  thermosonde  c  profiles  and  the  model. 
The  trends  are  similar  but,  of  course,  the  model  as  given  above  does  not  show 
turbulent  layering  structure,  and  it  cannot  be  expected  to  match  actual  turbulence, 
which  depends  on  the  Richardson  number  and  not  on  the  wind  speed  alone.  Further 

9 

studies  are  being  performed  to  compare  the  thermosonde  data  with  a  \  an  Zandt 
model,  which  does  depend  on  Richardson  number. 

4.4  Averaged  Thermosonde  C„^  Profile 

The  C2  data  for  the  10  day  and  9  night  flights  results  in  the  500  m  log-average 

profile  shown  in  Figures  51  and  52.  A  comparison  of  the  two  profiles  clearly  shows 

a  larger  C'2  in  the  daytime,  as  seen  in  Figure  53.  Starting  near  5  km,  the  difference 

increases  to  a  maximum  value  near  20  km  where  the  difference  is  four  times  the 
2 

nighttime  C'n.  Above  20  km,  the  difference  decreases,  finally  ending  above  34  km 
altitude.  A  careful  analysis  of  the  mierothermal  probes  (see  Appendix  A)  indicates 
the  diurnal  difference  cannot  be  attributed  to  solar  heating  of  the  probes. 

An  exponential  quadratic  model  may  be  used  to  fit  the  average  data.  The  model 
takes  the  form: 

log  C'2  =  A  +  Bz  +  Cz2. 


For  Z  2  3.5  km,  the  fitted  constants  are: 


1 0-Dav-flights 
Averaged 

A  -17.279 

H  v  o.  0852 

C  -  0. 0045 


9-  Night-flights 
A  ve  raged 

-16. 953 

-  0.  0251 

-  0.  0016 


9.  Van  Zandt,  T.  F.  ,  Gage,  K.  S.  ,  andtVarnock,  .J.M.  <  19811  An  Improved  Model 
for  the  Calculation  of  Profiles  of  C"  and  t  in  the  Free  Atmosphere  From 
Background  Profiles  of  W  ind,  Temperature  and  Humidity,  20th  Conference 
Radar  Meteor.  Soc.  .  Boston,  Amor.  Met.  Soc.  ,  pp  129-  135. 
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I'nr  summer  davtinic  conditions  at  the  experiment  coordinates,  the  model  sug- 

(|  •> 

e.-.t.-.  dropoff  rates  dli  km  10  -j—  log  C  )  of  -0.5  dii/km  and  -1.0  .IB- km.  at  a  km 

dz  ”  n 

nd  10  km  respeetivelv.  Nighttime  conditions  give  -0.4  dll  km  and  -1.2  dli  km  at 

9 

k.n  and  .'ill  km  res;  ti\rlv.  The  daytime  values  are  higher  than  the  night  values 

n  2 

nd  it  is  assumed  the  l  ister  daytime  decrease  of  C  with  altitude  would  eontinue 

n 

titil  trie  day  . nd  night  values  eoincide.  Thereafter,  they  would  both  den-ease  at  the 
ijh.'  rale.  If  it  is  assumed  that  t  .p  is  constant  with  altitude,  then  the  dli/ km  de¬ 
fer. so  e.m  he  interpreted  as  an  atmospheric  settle  height  (lit  using  the  relationship 


t>.  69 

11  r.TITTnTV 

I'he  observed  -1.  2  dli/km  tit  20  km  during  the  night  is  nearly  identical  to  the  standard 
atmosphere  pressure  scale  height  of  0.7  km  (-1.  2!>  dli/km!. 

4.5  Coherence  Scales 

2  1 
for  those  flights  where  thermosonde  C  exists,  the  transverse  coherence  length 

r  i  uttl  isopl  untie  angle  H  were  ealcuiated.  The  coherence  length  is 

I  .,  zf  .,  -2.5 

r(  1 .  o  j  k  /  ( ‘  (z)  dz  <m  1  . 

(  d 

and  the  isnplunatic  angle  ts 


f  (  ~  t/.t  dz 
1  n 


( rad)  . 


where  the  wave  number  k  is  211/ A,  z  is  the  observer  altitude  and  z.  is  the  maximum 

o  I 

altitude.  Results  of  these  calculations  are  shown  in  'fable  1.  Measurements  of  C“ 

n 


at  fl-12  m  intervals  were  logarithmically  interpolated  to  10  m  intervals  and  data  in 
the  noise  were  set  to  'noise'  =  0.  002’. 

For  those  flights  having  data  at  least  to  16  km,  the  total  log-average  of  the  iso- 
planatic  tingle  6 ^  is  calculated  for  day  and  night  conditions.  The  coherence  length  r^ 
could  be  calculated  including  the  two  low  altitude  flights.  The  averaged  values  are. 


r  =  6.7  _r 
o 

1.  2  cm 

llav 

TT  =  2.4: 
o 

0.  5  /i  rad 

r  =  7.6: 
o 

2.  4  cm 

Night 

0  -  5.  1  : 
o 

I.  I  u  rad 

where  the  errors  are  ±  1  standard  deviation.  As  expected,  daytime  conditions  pro¬ 
vide  smaller  values  of  the  coherence  length  and  thus  harsher  conditions  with  which 
a  pa rticula r  optical  system  must  contend. 


10.  fried,  L).  I..  (19061  ( iptical  resolution  through  a  randomly  inhomogeneous  medium 

for  -.erv  long  and  verv  short  exposures,  J,  Opt.  Soc.  Am.,  56:137  2-  1379. 

-  w'VV' 

11.  'A  .Iters,  11.1..,  Fnvier,  II.  I..,  and  Hines,  J.  R.  11979'  V  ertical  path  measure¬ 

ments  MTf  measurements,  .J.  ( )pt.  Soc.  Am.,  09:828-1137. 
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Figure  3.  Radiosonde  Ambient 
Temperature  and  Relative  Humidity 
Profiles  -  Flight  QX0009. 

Launched  at  17  30Z  on  21  May  1983. 
Humidity  data  are  truncated  above 
approximately  10  km 
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Figure  4.  Radiosonde  Ambient 
Temperature  and  Relative  Humidity 
Profiles  -  Flight  y\0050, 

Launched  at  2238/.  >n  30  dune  1983. 
Humidity  data  are  truncated  above 
approximately  10  km 
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Figure  5.  Radiosonde  Ambient 
Temperature  and  Helative  Humidity 
Profiles  -  F'light  QXOOO  1. 

Launched  at  0210Z  on  1  July  1983. 
Humidity  data  are  truncated  above 
approximately  10  km 
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Figure  6.  Radiosonde  Ambient 
Temperature  and  Relative  Humidity 
Profiles  -  Flight  QX0008. 

Launched  at  2033Z  on  11  July  1983. 
Humidity  data  are  truncated  above 
approximately  10  km 
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Figure  7.  Radiosonde  Ambient 
Temperature  and  Relative  Humidity 
Profiles  -  Flight  QX0045. 

Launched  at  0116Z  on  12  July  1983. 
Humidity  data  are  truncated  above 
approximately  10  km 
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Figure  8.  Radiosonde  Ambient 
Temperature  and  Relative  Humidity 
Profiles  -  Flight  QX0007. 

Launched  at  0228Z  on  16  July  1.983. 
Humiditv  data  are  truncated  above 
approximately  10  km 
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Figure  9.  Radiosonde  Ambient 
Temperature  and  Relative  Humidity 
Profiles  -  Flight  QX0036. 

Launched  at  0137Z  on  23  July  1983. 
Humidity  data  are  truncated  above 
approximately  10  km 


Figure  10.  Radiosonde  Ambient 
Temperature  and  Relative  Humidity 
Profiles  -  Flight  QX0012. 

Launched  at  2045Z.  on  23  July  1983. 
Humidity  data  are  truncated  above 
approximately  10  km 
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Figure  13.  Radiosonde  Ambient 
Temperature  and  Relative  Humidity 
Profiles  -  Plight  QX0035. 

Launched  at  2042Z  on  28  July  1983. 
Humiditv  data  are  truncated  above 
approximately  10  km 
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Figure  14.  Radiosonde  Ambient 
Temperature  and  Relative  Humiditv 
Profiles  -  Flight  QX0043. 

Launched  at  214  1Z  on  29  July  1983. 
Humidity  data  are  truncated  above 
approximately  10  km 
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Figure  17.  Radiosonde  Ambient 
Temperature  and  Relative  Humidity 
Profile  -  Flight  QX0020. 

Launched  at  0128Z  on  31  .July  1983. 
Humidity  data  are  truncated  above 
approximately  10  km 
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Figure  18.  Radiosonde  Ambient 
Temperature  and  Relative  Humidity 
Profiles  -  Flight  QX0040. 

Launched  at  033  1Z  on  7  Aug  1983. 
Humidity  data  are  truncated  above 
approximately  10  km 
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Figure  19.  Radiosonde  Ambient 
Temperature  and  Relative  Humidity 
Profiles  -  Flight  QX0031. 

Launched  at  0550Z  on  8  Aug  1983. 
Humidity  data  are  truncated  above 
approximately  10  km 
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Figure  20.  Winds  Aloft 
Profile  -  Flight  QX0024. 
Launched  at  1749Z  on 
15  May  1983 
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Figure  21.  Winds  Aloft 
Profile  -  Flight  QX0009. 
Launched  at  1730Z  on 
21  Alay  1983 
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Figure  29.  Winds  Aloft 
Profile  -  Flight  QX0038. 
Launched  at  2044Z  on 
25  July  1983 
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Figure  30.  Winds  Aloft 
Profile  -  Flight  QX0035. 
Launched  at  2042Z  on 
28  July  1983 
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Figure  33.  Winds  Aloft 
Profile  -  Flight  QX0034. 
Launched  at  2016Z  on 
30  July  1983 
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Figure  34.  Winds  Aloft 
Profile  -  Flight  QX0020. 
Launched  at  0I28Z  on 
3  1  July  1983 
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Table  1.  Balloon  Flight  Parameters 


Launch 

Date 

(UT> 

Time 

Flight 

ID 

(rrt2  sec"2) 

r 

o 

(cm) 

(p  rad) 

Max. 

Alt. 

(km) 

Max. 

Wind 

Alt 

(km) 

Day  / 
Night 

23  Mar 

1850 

QX21 

- 

6.7 

2.6 

27.6 

- 

D 

15  May 

1800 

QX24 

389. 

5.  r 

2.  2 

29.  6 

26.  6 

D 

21  May 

1730 

QX09 

- 

- 

- 

32.  9 

32.4 

D 

30  Jun 

2237 

QX50 

163. 

8.  2 

3.  1 

24.  9 

25.  8 

D 

1  Jul 

0210 

QX01 

- 

- 

- 

27.  6 

- 

N 

11  Jul 

2033 

QX08 

- 

- 

- 

3.  4 

4.  0 

D 

12  Jul 

0116 

QX45 

1250. 

5.4 

3.  2 

19.  5 

8.  0 

N 

16  Jul 

0229 

QX07 

- 

- 

- 

21.0 

22.  6 

N 

23  Jul 

0137 

QX36 

1249. 

10.  3 

6.  7 

28.  4 

28.  7 

N 

23  Jul 

2045 

QX  1 2 

415. 

7.  6 

2.  6 

31.9 

32.  5 

D 

24  Jul 

0236 

QX  1 9 

- 

- 

- 

19.  1 

16.  5 

N 

25  Jul 

2044 

QX38 

484. 

5.4 

2.  1 

34.  5 

33.  1 

D 

28  Jul 

2042 

QX35 

619. 

7.4 

- 

7.  4 

7.  3 

D 

29  Jul 

2141 

QX43 

1805.  : 

7.  5 

- 

10.  1 

10.  3 

D 

30  Jul 

0113 

QX  1 1 

1256. 

7.  8 

3.  8 

24.  4 

23.  9 

N 

30  Jul 

2016 

QX34 

472. 

5.0 

1.6 

30.  0 

31.  1 

D 

31  Jul 

0128 

QX20 

469. 

10.  7 

5.  7 

31.  0 

31.0 

N 

7  Aug 

0331 

QX40 

179. 

5.  5 

4.5 

22.  9 

22.  2 

N 

8  Aug 

0549 

QX31 

158. 

6.  1 

4.  7 

23.  8 

25.  8 

N 

mean  square  wind  speed  from  5  to  20  km  (or  maximum  altitude  ) 
coherence  length 
isoplanatic  angle 
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Figure  37.  Log -Average  C 
Thermosonde  Profiles  and  n  the 
Hufnagel  Model  -  Flight  QX0021. 
Launched  at  18502  on  23  Mar  1983. 
Thermosonde  values  were 
averaged  every  500  m,  and  curves 
showing  the  mean  multiplied  by, 
and  divided  by,  1  sigma  values 
are  shown  on  either  side  of  the 
mean 
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Figure  38.  Log- Average  C 
Thermosonde  Profiles  and  n  the 
Hufnagel  Model  -  Flight  QX0024. 
Launched  at  1749Z  on  15  May  1983. 
Thermosonde  values  were 
averaged  every  500  m,  and  curves 
showing  the  mean  multiplied  by, 
and  divided  by,  1  sigma  values 
are  shown  on  either  side  of  the  2 
mean.  The  thick  line  shows  the  C 
values  from  the  Hufnagel  model  n 
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Figure  39.  Log-Average  C 
Thermosonde  Profiles  and  nthe 
Hufnagel  Model  -  f'light  QX0050. 
Launched  at  2238Z  on  30  June  1983. 
Thermosonde  values  were 
averaged  every  500  m,  and  curves 
showing  the  mean  multiplied  by, 
and  divided  by,  1  sigma  values 
are  shown  on  either  side  of  the  „ 


The  thick  line  shows  the  C 


values  from  the  Hufnagel  model 
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Figure  40.  Log-Average  (j 
Thermosonde  Profiles  and  nthe 
Hufnagel  Model  -  Flight  QX0045. 
Launched  at  0116Z  on  12  July  1983. 
Thermonsonde  values  were 
averaged  every  500  m,  and  curves 
showing  the  mean  multiplied  by, 
and  divided  by,  1  sigma  values 
are  shown  on  either  side  of  the  „ 


The  thick  line  shows  the  C 


values  from  the  Hufnagel  model 
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Figure  41.  Log-Average  C 
Thermosonde  Profiles  and  n  the 
Hufnagel  Model  -  Flight  QX0036. 
Launched  at  0137Z  on  23  July  1983. 
Thermosonde  values  were 
averaged  every  500  m,  and  curves 
showing  the  mean  multiplied  by, 
and  divided  by,  1  sigma  values 
are  shown  on  either  side  of  the  2 
mean.  The  thick  line  shows  the  C 
values  from  the  Hufnagel  model  n 


2 

Figure  42.  Log-Average  C 
Thermosonde  Profiles  and  n  the 
Hufnagel  Model  -  Flight  QX0012. 
Launched  at  2045Z  on  23  July  1983. 
Thermosonde  values  were 
averaged  every  500  m,  and  curves 
showing  the  mean  multiplied  by, 
and  divided  by,  1  sigma  values 
are  shown  on  either  side  of  the  2 
mean.  The  thick  line  shows  the  C 
values  from  the  Hufnagel  model  n 
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Figure  43.  Log-Average  C 
Thermosonde  Profiles  and  1  the 
Hufnagel  Model  -  Flight  QX0038. 
Launched  at  2044Z  on  25  July  1983. 
Thermosonde  values  were 
averaged  every  500  m,  and  curves 
showing  the  mean  multiplied  by, 
and  divided  by,  1  sigma  values 
are  shown  on  either  side  of  the  ^ 
mean.  The  thick  line  shows  the  C 
values  from  the  Hufnagel  model  n 


Figure  44.  Log-Average  C 
Thermosonde  Profiles  and  n  the 
Hufnagel  Model  -  Flight  QX0035. 
Launched  at  2042Z  on  28  July  1983. 
Thermosonde  values  were 
averaged  every  500  m,  and  curves 
showing  the  mean  multiplied  by, 
and  divided  by,  1  sigma  values 
are  shown  on  either  side  of  the  ^ 
mean.  The  thick  line  shows  the  C 
values  from  the  Hufnagel  model  n 
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Figure  47.  Log-Average  C 
Thermosonde  Profiles  and  n  the 
Hufnagel  Model  -  Flight  QX0034. 
Launched  at  2016Z  on  30  July  1983. 
Thermosonde  values  were 
averaged  every  500  m,  and  curves 
showing  the  mean  multiplied  by, 
and  divided  by,  1  sigma  values 
are  shown  on  either  side  of  the  2 
mean.  The  thick  line  shows  the  C 
values  from  the  Hufnagel  model  n 


Figure  48.  Log-Average  C 
Thermosonde  Profiles  and  n  the 
Hufnagel  Model  -  Flight  QX0020. 
Launched  at  0128Z  on  31  July  1983. 
Thermosonde  values  were 
averaged  every  500  m,  and  curves 
showing  the  mean  multiplied  by, 
and  divided  by,  1  sigma  values 
are  shown  on  either  side  of  the  2 
mean.  The  thick  line  shows  the  C 
values  from  the  Hufnagel  model 
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Figure  49.  Log-Average  C 
Thermosonde  Profiles  and  n  the 
Hufnagel  Model  -  Flight  QX0040. 
Launched  at  033 1Z  on  7  Aug  1983. 
Thermosonde  values  were 
averaged  every  500  m,  and  curves 
showing  the  mean  multiplied  by, 
and  divided  by,  1  sigma  values 
are  shown  on  either  side  of  the 
mean.  The  thick  line  shows  the  C 
values  from  the  Hufnagel  model 
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Figure  50.  Log-Average  C 
Thermosonde  Profiles  and  n  the 
Hufnagel  Model  -  Flight  QX0031. 
Launched  at  0550Z  on  8  Aug  1983. 
Thermosonde  values  were 
averaged  every  500  m,  and  curves 
showing  the  mean  multiplied  by, 
and  divided  by,  1  sigma  values 
are  shown  on  either  side  of  the 
mean.  The  thick  line  shows  the  C‘ 
values  from  the  Hufnagel  model  1 
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Figure  51.  Log-Average  C 
of  Ten  Daytime  ThermosoncPe 
Flights 
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Figure  52.  Log-Average  Cn 
of  Nine  Nighttime  Thermosonde 
Flights 
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Figure  53.  Normalized 
Difference  Between  Day 
and  Night  Measurements  of 
r2  =  , Pay- Night  i 
^n  1  Night  ’ 
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5.  THERMOSONDE  -  RADAR  COMPARISON 


2  2 

For  those  flights  where  radar  data  exists,  composite  radar  Cn  profiles  were 

constructed  by  assembling  the  data  in  500-m  bins.  It  was  desired  to  compare  the 

radar  data  with  the  500-m  log  averaged  thermosonde  data  as  closely  as  possible  in 

time  and  space.  Since  the  log-average  thermosonde  data  is  centered  at  known  times 

and  altitudes,  radar  data  could  be  placed  into  altitude  bins  corresponding  to  the 

balloon  time  ±  250  sec  and  altitude  ±  250  m.  This  being  done,  the  data  in  each  bin 

were  log-averaged  and  resulted  in  those  profiles  shown  in  Figures  54-61.  Since  each 

individual  raw  radar  profile  was  provided  for  each  minute  of  flight,  the  points  on 

2 

the  graphs  are  log-averages  of  about  9  radar  values.  However,  due  to  changes 
in  the  balloon  ascent  rate  and  other  factors,  the  actual  number  of  points  varies 
between  5  and  12. 

Considering  the  difference  in  resolution  between  the  two  data  sets  and  the  re¬ 
maining  separation  in  horizontal  position  between  the  measurements,  the  agreement 
is  good.  Figure  62  shows  plots  of  the  mean  differences  between  the  radar  and 

thermosonde  turbulence  for  six  combined  flights.  Generally,  on  average,  the  radar 

2  2 
measures  larger  C’n  by  about  1.  2  to  2.  0  times  the  thermosonde  C'n. 
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Figure  54.  Comparison  of 
500-m  Log-Average 
Thermosonde  Profiles  and 
Composite  Log-Average 
Radar  Profiles  -  Flight  QX0050, 
Launched  at  2238Z  on 
30  June  1983 
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Figure  55.  Comparison  of 
500-m  Log-Average 
Thermosonde  Profiles  and 
Composite  Log-Average 
Radar  Profiles  -  Flight  QX0045. 
Launched  at  0116Z  on 
12  July  1983 
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Figure  56.  Comparison  of 
500-m  Log -Average 
Thermosonde  Profiles  and 
Composite  Log-Average 
Radar  Profiles  -  Flight  QX0012. 
Launched  at  2045 Z  on 
23  July  1983 


Figure  57.  Comparison  of 
500-m  Log-Average 
Thermosonde  Profiles  and 
Composite  Log-Average 
Radar  Profiles  -  Flight  QX0035. 
Launched  at  204  2Z.  on 
28  July  1983 
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Figure  60.  Comparison  of 
500-m  Log-Average 
Thermosonde  Profiles  and 
Composite  Log-Average 
Radar  Profiles  -  Flight  QX0034 
Launched  at  2016Z  on 
30  July  1983 
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Figure  61.  Comparison  of 
500-m  Log-Average 
Thermosonde  Profiles  and 
Composite  Log-Average 
Radar  Profiles  -  Flight  QX0020. 
Launched  at  0128Z  on 
31  July  1983 
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Figure  62.  Plot  of  Mean 
Differences  Between  Thermosonde 
and  Radar 
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6.  THERMOSONDE  -  RADAR  -  SCINTILLOMETER  COMPARISON 

On  the  evening  of  30  July  1981,  the  thermosonde  (QX0020),  radar,  and  scin- 

2 

tillometer  operated  together  successfully.  Results  of  their  C  measurements  are 
shown  in  Figure  63.  The  thermosonde  profile,  shown  by  the  thin  line,  starts 
0.  5  km  above  the  surface  and  ends  at  30  km.  This  is  a  500-m  weighted  log-average 
of  the  data.  The  radar  profile,  shown  by  asterisks,  is  a  time/altitude  weighted 
log-average  composite  of  the  data.  It  tends  to  match  radar  time  and  altitude  with 
the  balloon  trajectory,  though  not  in  the  horizontal  plane.  The  scintillometer 
profile,  shown  by  the  thick  line,  was  observed  at  0208  h  UT  when  the  balloon  was 

at  12.5  km.  Over  a  1-h  time  period,  the  maximum  excursion  of  the  scintillometer 

2  3 

C  data  at  any  one  altitude  was  approximately  1  order  of  magnitude.  The  trend 

line  of  the  thermosonde  profile  coarsely  agrees  with  this  one  scintillometer  profile, 

but  until  additional  comparisons  are  made  only  the  suggestion  of  agreement  is 

postulated.  More  extensive  thermosonde/scintillometer  comparisons  were  conducted 

in  March  1983  in  Boulder,  Colorado.  These  will  be  reported  in  a  future  report. 
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Figure  63.  Comparisons  Between  Log-Average 
Thermosonde.  Radar,  and  Scintillometer  cj" 
Profiles  -  One  Flight 


7.  THERMOSONDE  -  AIRCRAFT  COMPARISON 

2 

Hot  wire  anemometer  experiments  were  conducted  from  the  AKA  airplane 
on  11  July  1981  and  16  July  1981.  Thermosonde  and  airplane  experiments  were  con¬ 
ducted  together  successfully  on  the  night  of  11  July.  The  thermosonde  launch  was  at 
0115  h  UT  while  the  aircraft  took  measurements  from  0133  to  0247  h  LT.  Measure¬ 
ments  were  taken  during  ascent  to  8.  5  km  and  then  during  descent  to  0.  5  km.  The 

resulting  profiles  are  compared  to  the  thermosonde  measurements  in  Figures  64-67. 

2 

Two  separate  bridge  circuits  measured  aircraft  C  data  independently.  Bridge  No.  1 
2  ^ 

measured  C  values  about  2  to  4  times  that  of  the  thermosonde  except  for  the  lowest 

altitude  descent  data.  Bridge  No.  2  was  high  by  an  order  of  magnitude  and  was  not 

correctly  calibrated.  Further  evidence  of  aircraft  measurement  error  was  seen  in 

comparative  measurements  of  27  Dec  1983  and  30  Dec  1983  where  the  aircraft  showed 

2 

1  to  2  orders  of  magnitude  higher  C  .  Clearly,  the  ARA  aircraft  measurements  pro¬ 
duce  gross  error.  This  technique  is  not  recommended  for  future  programs  unless 
more  study  produces  better  information  on  the  sources  of  the  error. 
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Appendix  A 

Solar  Heating  Effects  on  Probe  Wire 


To  determine  the  extent  of  possible  thermosonde  measurement  error  due  to 
solar  heating  of  the  probe  wires,  the  following  heat  transfer  analysis  is  suggested. 
External  solar  heating  of  the  wire  is  seen  to  occur  by  direct  solar  radiative  heating. 
Egfy  by  solar  radiative  reflection  of  the  atmosphere,  Eg^,  and  by  long  wave  radia¬ 
tion  from  the  earth,  E^.  Probe  cooling  is  seen  to  occur  by  long  wave  radiation 
from  the  wire  E^y  and  by  convective  cooling,  E^. 

Estimates  of  the  heat  transport  are  obtained  from  the  following: 


fS  t  _  A 
r  s 


"SD 


where 


f  =  estimate  of  fraction  of  soLar  radiation  reaching  a  given  altitude 
f  *  0.  5  at  surface 
f  =  0. 9  at  30  km 

-2 

S  =  solar  constant  =  0.  14  W  cm 

A  =  surface  area  of  wire  2  Jt  RL 
s 

€ rj.  estimate  of  wire  absorbitivity  for  visible  wavelength  ~0.4 
F._  -  form  factor  for  direct  radiation  =  n 


where 


a  =  estimate  of  fraction  of  reflected  radiation  =  0.  35 
=  form  factor  for  reflected  radiation  =  2.  0 

S  e  T  ...  A 
e  L\V  s 


where 


S  =  estimate  of  heat  flux  from  earth  at  a  given  altitude 


S  ~  0.  066  W  cm  at  surface 
e 

S  =  0.  022  W  cm  2  at  30  km 
e 


e  =  estimate  of  wire  emissivity  for  long  wave  radiation  =  0.  03 


ELW  =  SLWAS  aTs  StLW  As  °(Tg  +  4  AT  ] 


where 


a  =  Stephan- Boltzman  constant  =  5.67  X  10  12  W  cm-2  K~4 


Tg  =  wire  temperature,  K 

Tg  =  air  temperature,  K  =  300  K  at  surface  and  225  K  at  30  km 
AT  =  (Tg4  -  Te) 

T4  *  IT  +  4  AT  T3  for  AT  <  <  T 
s  e  e  e 


E  =  J  A  hAT 
c  s 


where 


=  mechanical  equivalent  of  heat  =  4.  18  W/cal/s 
=  estimate  of  convective  heat  transfer  coefficient 

-2  - 1  _ j 

*  0.  12  cal  cm  sec  K  at  surface 
a  0.  03 1  cal  cm  2  sec  *  K  *  at  30  km . 


Now  since,  Eg^  +  Eg^  +  E^^,  =  E^yy  +  Ec  anc*  combining  terms. 


for  daytime  conditions 


AT  =  0.  03  deg  at  surface 
=  0.  19  deg  at  30  km 

for  nighttime  conditions,  the  solar  radiative  term  vanishes,  and 


AT  =  0  . 


The  conclusions  are  of  course  that  we  need  not  be  concerned  with  nighttime 
heating.  Daytime  solar  heating  will  not  raise  the  wire  temperature  more  than 
0.  03  deg  and  0.  2  deg  above  ambient  conditions  at  the  surface  and  30  km  respectively. 
Since  this  is  a  dc  effect  upon  both  wires,  the  ac  bridge  will  remove  the  gross  im¬ 
balance,  if  any.  Small  velocity  fluctuations  in  the  air  flow  over  the  wires  can  cause 
small  variations  in  the  heat  transfer  coefficient,  especially  at  the  surface.  A  large 
fluctuation  in  windspeed,  say  10  percent,  can  cause  at  most  a  10  percent  change  in 
the  heat  transfer  coefficient.  This  will  produce,  at  most,  an  ac  temperature  fluc¬ 
tuation  of  0.  003°  C,  an  undetectable  fluctuation.  The  velocity  effect  on  the  heat 
transfer  coefficient  falls  off  rapidly  with  altitude.  At  10  km,  a  10  percent  velocity 
fluctuation  can  cause  very  little  temperature  fluctuation.  At  30  km,  no  change  in  the 
heat  transfer  coefficient  will  be  caused  by  forced  convection  and  hence  no  detectable 
temperature  fluctuation  will  be  produced. 


